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OBSERVATIONS OF INFALL AND CAUSTICS 

R. Brent Tully 1 



Abstract. Identification of first and second turnaround radii in groups 
of galaxies provide interesting scaling relations and may constrain cos- 
mological parameters. 



1 Introduction 
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, Collapse timescales, t c , depend on the density, p, of a bound structure which in 
C^S ' turn depends on the mass, M, enclosed within a radius, r: t c ~ p~ x / 2 ~ (r 3 /M) 1 / 2 . 
, For two structures at the same phase of collapse today: 

n = r a (M 1 /M 2 ) 1/3 . (1.1) 

O !' Suppose we are comparing structures within the radius of second turnaround, 
q . T^ti the- maximum radius that a collisionless object reaches upon falling once 
5^ ■ through the center of the group potential. The dynamical stages of the two com- 
"ry^ \ parison groups are similar so M ~ o-yr 2 t where cry is the velocity dispersion 
& ' within the quasi- virialized region bounded by r%t- If we consider t = today, then 
>• (< J v r ' 2t l r 2t) 1 ^ 2 = constant and: 

' ' " ay^r 2t . (1.2) 

^: 

>r 2 Observations of the caustic of second turnaround 

Models of spherical infall of cold, collisionless particles (Bertschinger 1985, ApJS, 
58, 39) suggest that the radius of second turnaround r 2 t will be marked by a density 
drop and a transition from a large velocity dispersion interior to r 2 t to cold infall 
exterior to r 2t . Projection effects in the real world might make the cold infall 
difficult to distinguish, but that is a matter to be explored below. The theoretical 
models also predict density cusps at the 2nd and all interior turnarounds but it 
can be presumed that these features would be wiped out in the real world with 
asphericity, finite peculiar motions, and collisions. 
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It is not surprising that the caustic of 2nd turnaround has received little obser- 
vational attention because the anticipated density transition is subtle and velocity 
effects are obscured by projection. Fortunately, in a recent study (Mahdavi et al. 
2005, AJ, astro-ph/0506737) a particularly clean case came to light. The modest 
but dense group of early-type galaxies around NGC 5846 lies within the Local 
Supercluster at 26 Mpc in relative isolation. The group is a knot within a filament 
that lies close to the plane of the sky and the space to the foreground and back are 
voids. The group has a huge dwarf population with 325 candidate members and, 
overall, velocities are available for 85 affirmed members. The NGC 5846 group has 
a velocity dispersion of 322 km s _1 and a mass of 8 x 10 13 M©. 
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Fig. 1. Radial gradients. Top: Surface density dependence on radius for galaxies with 
observed velocities and R < 17.2. The surface density drops smoothly with radius out 
to ~ 1.8° then drops abruptly to a roughly constant value. Bottom: Velocity dispersion 
as a function of radius. Each data point represents the velocity dispersion of all galaxies 
within 0.5 Mpc of individual galaxies considered one at a time; thus, the data points are 
not independent. The sloping line within 1.8° is fit to velocity dispersions averaged over 
the 4 annular bins of the top panel. The flat line, showing the mean velocity dispersion 
outside 1.8°, lies at 66 km s~ . 

Figure 1 summarizes the pertinent observations. Sloan Digital Sky Survey 
velocities are available for all galaxies in the region with R < 17.2 which provides 
a complete magnitude limited sample uncontaminated by background. The top 
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panel of Fig. 1 shows the run of surface densities and the bottom panel shows the 
run of velocity dispersions. The latter, to be precise, is the velocity dispersion in 
small windows centered at each galaxy in the sample. The point is to show the 
very small dispersions in local regions at large radii, the signature of cold infall. 
If dispersions were averaged over annuli then the cold flow at large radii would 
be camouflaged by projection differences between widely separated galaxies. The 
anticipated transition at r<u is seen at 1.84° = 840 kpc. 
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Fig. 2. Correlation between velocity dispersion and inferred radius of 2nd turnaround 
for E/S0 knots in Local Supercluster. 

The locations of 2nd turnaround are less certainly defined for other nearby 
knots of E/S0 galaxies but rough estimates are available. Those estimates give 
rise to the correlation seen in Figure 2. The straight line describes the relation 
anticipated by Eq. 1.2: 

cry /r 2 t = 390 km s" 1 Mpc" 1 (2.1) 
3 The ratio of first to second turnaround radii 

The ratio of the radii of first and second turnaround, Txtl^iu depends on cos- 
mology. Mamon (private communication) finds the related ratio r u /r V i r ranges 
between 3.1 and 3.7 for < fl m < 1 in a flat Universe. The best shot at defining 
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r\t with present data lies with the Local Group because in this unique circum- 
stance we see almost the entire vector of radial infall motion. The radii r 2 t around 
the Milky Way (MW) and M31 can be inferred from Eq. 2.1 or from the vari- 
ant that follows from Eq. 1.1. Given mass estimates for these galactic systems of 
1 x 10 12 M Q (Wilkinson & Evans 1999, MNRAS, 310, 645; Evans et al. 2000, ApJ, 
540, L9) then r^t ~ 200 kpc. It is seen in Figure 3 that 200 kpc is roughly the 
domain of the gas poor and high velocity dispersion galaxies around M31 and MW. 
Beyond this radius, galaxies are mostly gas rich and have velocities indicative of 
infall. The transition to blueshifts, marking the surface of first turnaround, occurs 
at a radius r\ t ~ 900 kpc. In principle, the ratio ri t /r2t provides a constraint on 
Q m . The present observed parameters are soft because of the small number of 
Local Group galaxies and the dumbbell shape of the potential. In any event, it 
would be nice to see what simulations anticipate in similar circumstances. 
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Fig. 3. Discrete quasi-virialized and infall regions within the Local Group. Velocities 
with respect to the centroid of the Local Group are plotted against the distance of a 
galaxy from the nearer of M31 or the Milky Way. Locations of 2nd turnaround for M31 
and M31 are scaled using Eq. 1.1 from the relation for E/S0 knots shown in Fig. 2. The 
approximate position of the 1st turnaround for the combined M31+MW system is also 
indicated. 



